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ABSTRACT 


DISPERSION OF DYES AND POLLUTANTS IN NORTH 


SASKATCHEWAN RIVER 


The object at ehussthesis .wasetoustudvyethe Capacity 
of North Saskatchewan river to assimilate wastes, parti- 
aid fem: the problem of dispersion of pollutants. The study 
area involved the reach of river below Edmonton to the 
Vinceasbridge. 

Field measurements and sampling were carried out 
during the winter months from September Sik Seren MER acelele 
W972, becauceuthisets thesperiod Of tow flows io) LReeri ver 
and oe ice cover. The examination of the river water took 
place regularly at two locations ~- Fort Saskatchewan bridge 
and Vinca bridge sampling stations. At both locations sev — 
eral samples were taken across the river cross-section to 
determine the distribution of pollutants. The results of 
chemical laboratory tests confirmed the suspected incomplete 
mixing below the Town of Fort Saskatchewan, whereas the Fort 
Saskatchewan bridge cross-section showed quite uniform 
distribution Of pollutants. 

In January 1972 an extensive survey of eletey Ve Eiete 


cross-section at Vinca was carried out in order to obtain 
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iv 
the necessary data for a dispersion coefficient: determination. 
This, together with a dye tracer study, provided data for 
mathematical solution of the dispersion equation, applic- 
ap lemEOnril ce Teaci@eOreviVver. 

A twenty nareene solution of Rhodamine W. T. was 
used in the dye tracer study. Two runs were made, the 
first, from Dow Chemical at Fort Saskatchewan to Vinca, 
showed incomplete mixing at Vinca and provided the basis 
for the development of the three-dimensional model of dis- 
persion, the second, from Edmonton to Fort Saskatchewan, 
showed quite uniform mixing at Fort Saskatchewan Bridge. 

The results of the three-dimensional model of dis- 
pension as well as the first dye test was the basis for 
determination of coefficients relating the mean concentra- 
tion in the river to the values obtained at the Government 
of Alberta, Department of Environment sampling point. 

These coefficients, shown in the last chapter, indicate the 
dispersion in the stretch of river between Fort Saskatchewan 


and Vinca. 
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GLOSSARY OF TERMS 


ASSIMILATIVE CAPACITY 


DISPERSION 


LEFT BANK (RIGHT BANK) 


WATER POLLUTION 


The capacity of a river to assimi-— 
late waste loads. 

The mixing and spreading of con- 
centrated mass of pollutant or dye 
tracer in a Natural) stream. 
Looking downstream the right and 
left bank is considered to be the 
Onewon the right andslertt side 
subsequently. 

Impairment of water quality and 
the resulting adverse effect on 


water use. 
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CHAP Rar 
INTRODUCTION 


ieee eOUlLUtI ONsOreNaturaleotreams. 

Not the least of the important uses of water is its 
use for the disposal of liquid wastes from cities and in- 
dustries. There is nMothing on. the horizon to suggest the 
discovery of waste handling systems superior in efficiency 
and economy to our present water carriage systems. We 
will continue to discharge wastewater into our streams. We 
havemnomcnolccsinethisemacter slum only choice lies in the 
extent to which we allow the discharge of wastewater of all 
saves to impair and interfere with other beneficial water 
uses. To the extent that we can eliminate or redauce water 
quality damage and consequent interference with water use, 


we have accomplished water pollution control. 


2. Assimilative Capacity of) Natural Streams. 

With ever increasing demands being placed on surface 
waters for disposal of man's waste products, it is essential 
that the capacity of these waters for assimilating wastes be 
known for the protection of their many beneficial uses. The 
determination of this assimilative capacity is dependent on 


three important factors: 
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1. The dispersion of the wastes and the dilution 
effected; 
2. The net advection of the waste receiving water 
system; 
3. The decomposition of the wastes. 
To permit the prediction of the concentrations of waste 
Bede itents in receiving waters, all three factors either 
must be known or estimated realistically. The natural 
ability of rivers to assimilate waste is limited and varies 
from stream to stream. This variation is largely dependent 
on character of wastes, eral cClamaticecondveions, amount 


of waste and character of the stream. 


oe Past Probilemseand Research on North Saskatchewan River. 
The examination OL the river water in the early 
1950's showed very low levels of dissolved oxygen down- 
stream of Edmonton. In 1955-1956 it was found that the 
oxygen in the river at Duvernay, Lindbergh and Lloydminster 
was depleted and septic conditions which developed in ab- 
sence of oxygen produced strong odors. This situation re- 
sulted in an attempt by the Alberta Department of Public 
Health to artificially aerate the river. The reaeration 
was unrortinately, not Ssuccessiul., “in the rollowing years 


more stringent controls of all waste discharges were 
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3 
introduced and oxygen levels in the river during the winter 
months were significantly higher. The situation was fur- 
ther improved by the new sewage treatment plant at Edmonton, 
completed in 1956 for primary treatment and 1957 for sec- 
ondary. Newly ore erent ed sewage lagoons for packing house 
wastes in 1965 also reduced the load of wastes in the winter 
months. 

Hie \oCasuMmmMeCtOLeloy mine ebiterndisciplinary. Commit— 
tee for Environmental Quality (ICEQ) and the Fish and Game 
Association, together with the Edmonton Anti-Pollution 
Greup, conducted@a sunvey and published jointly the study 
"Effluent Flow Patterns in the North Saskatchewan River" 
Go). “WDE study attempted to determine the extent of 
mixing of effluents discharged into the North Saskatchewan 
River in the Edmonton-Fort Saskatchewan area. It also 
pointed out the fact, that sampling on the river is usually 
done by taking 'grab' samples at midstream, whereas almost 
all effluents in the Edmonton area enter the river on the 
south bank. The Government of Alberta, Department of 
Environment monitor, installed every year at Vinca Bridge, 
has a pipe intake extending about 150 feet from the north 
bank and is, according to the Edmonton Anti-Pollution Group 


study, providing samples which do not reflect the true 
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river water quality. 


4. Study Objectives. 

Most studies of the natural purification in streams 
have been carried out on rivers of more southerly lati- 
tudes, flowing under warm conditions. The Department of 
Civil Engineering, at the University of Alberta has been 
involved for several years in studying and investigating 
pollution problems, biooxidation and oxygen utilization in 
the North Saskatchewan River under ice cover at tempera- 
CuGnesecloses tou degrees Ce) This study is a further step 
tnechnas continuingsprogramaand 1ts object is to study 
dispersion in the North Saskatchewan River under ice cover 
so as to aid in the selection of representative sampling 
points. 

DiC GOdLsOLMULniS study. is >: 

1) To determine a dispersion coefficient for the 
dispersion equation, which may be applied to 
pollutants entering the river at or above Fort 
Saskatchewan and measured at the Vinca Bridge 
Government sampling station. 

2) To determine the concentration distribution of 
wastes at Fort Saskatchewan Bridge. 


3) To determine the concentration distribution of 
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effluents from Fort Saskatchewan in the river 
cross-section at Vinca. 

To determine the relation between the. mean concen- 
tration in Vinca cross-section and the values 
obtained at the point of government sampling. 

To verify the distribution by means of dye 


tracer test. 
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CHAPTER II 
STUDY AREA 


1. General Discription. 

The North Saskatchewan River = parteor the Hudson 
Bay drainage system, rises in the Rocky Mountains near the 
continental divide. The high summer flow is derived from 
the melting snow and glaciers in the mountains, whereas the 
low winter flow is mainly from lakes and groundwater stor- 
age. 

The river flows in an easterly direction through 
Alberta to Saskatchewan, where it is joined by the South 
Saskatchewan River, “From this junction it flows again in 
an easterly direction to Manitoba and into Lake Winnipeg. 
The upper part (upstream of Drayton Valley) drops from an 
elevation about 5000 feet to 2400 feet in about 200 miles. 
Below Drayton Valley the slope diminishes so that at 
Edmontcn the Water elevation is 2010 feet above the geodetic 
datum and 210 miles downstream at the Alberta-Saskatchewan 
border it drops to an elevation of 1639 feet. in our 
study area the actual drop is 72.6 feet between Edmonton 
aHaGaVinca ferry, a distance Gr 39.8 miles. 


There are only small tributaries in this stretch of 
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river and they are almost completely frozen during the 
winter months; therefore, the rate of flow is uniform be- 
tween Edmonton and Vinca ferry. The river is covered with 
PCeRECOCMO CO OemMOncNs cacheycar. 

Tce formation Ustally appears by the end of October 
or beginning of November and remains until April. Duration 
of ice cover in the period 1915-1965 and its frequency 
distribution was studied by R. D. Cameron (1) and Figure l 


mo takens trom his Study 


2. River Dascharge. 

The mountain and foothill portion of the North Saska- 
tchewan River basin provides a greater part of the stream- 
flow than do the plains’ of the prairie provincies,. The 
average yield of the higher back range areas is over 40 
inches per year, while the drier plains yield less than one 
inch. In relatively dry years, such as 1948-1949, the 
plains portion of the basin with 87% of the area Conti. 
buted less than 16% of the flow (see Table I). 

There are large variations in the annual maximum and 
minimum daily discharge. There is a tendency toward two 
peaks a year, but not necessarily in every year. The fLcst 
maximum occurs as a result of snowmelt and rainfall in 


April and May, the second, in June or July, is due to 
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FIGURE 2 : Map Showing the Study Area. 


(Arrows indicate sampling points and points of release of dye) 


& by, 
peer 


7) 


ae o & 
a Re oe 





¢ 





‘ 


10 





GST 0°OOT 000’0O0T6‘E 0°O0OT 00S‘LY TWLOD 
GAO LarG:h 000’9Z29 €°L8 SLU OLY SUuTeTd 
Tac Cunt 000‘07 G0 aA STTty}OOg TeyIO 
ZSeAOT ASZEMACSTO 
T°8 Lesy 000°rPT6‘T v°6 890‘ SATOSOY WSerTOg 
0°02 eer 000‘0S6 6°T T68 yjueg 
Ceael Sh 2s: 000‘0SE 6°0O SPP aZedsep 
aot G7 S 000‘00€/T spore ores 9 Syreq TPUOTIEN 
soyoul MOTIWeCSIAS qZeSe75 o20V eolay scp DS 


eoazy eomnos 
PTETA TOEREM TROL FO ¢ MoTyweet4zS uTSeg FO 3 erolyv v = 


A 


676T-876T - YHATY NVYMAHOLVASVS HLAYON GHG AO MOTANWEULS 


I dIdvi 








7 rr ~ - = = 7 
A. Woltgeved2 wigtt Sp? 
a Jaga ston. , "eam 


———— eS — 2—————— 
all . 

ris “ ~ ” r 
Acie *“ i ok be 


prt - Bio O0G ce e.i) : F 


} 
pe 


- Scuet. - ests oTeaet ¢.2 
e : . 
7 4 sO . aa i 
78: >: i. WORE. be Bead 
= 
Blt 47 G00 1.0 P rex 
“ 
s mel | 7 
r { Othe r ~~ | ee: 
:& 
‘Gl ‘a r 
c U sie 
agi ———— a ~ —_— — _ — —— — ee — — —~ 


of 


12 


MAIN EFFLUENT DISCHARGERS INTO NORTH SASKATCHEWAN 


RIVER IN ALBERTA 








Character Type 
Source of of 
Waste Treatment 
Town of Rocky Mountain House Domestic none 
Town of Drayton Valley Domestic lagoon 
Town of Devon Domestic ST 
Imperial Oil - Devon Industrial Lagoon 
City of Edmonton Domestic ST in winter 
Sedan Industries Ltd. Industrial well 
Building Products Ltd. Industrial AL 
Tuperial Omeutdr industria iN DAIL Ag, Wee 
Texaco tindus tie. al Papel, weLl 
Union Carbide Industrial lagoon 
Ses. ls Oils Industrial lagoon 
Gulf Oil Canada Ltd. Industrial Meith ibe, 7, aWeleleleys 
Chemcell Industrial well, lagcon 
Uniroyal Industrial lagoon 


Alberta Hospital - Oliver Waste stored in lagoons 
City of Edmonton Packing Plants during winter months and 


Sherwood Park discharged in summer 


cr EN A Sm 


si. 


WAWSNITSA Ss 


ee 
. oe] 
al 
uf \ + 
= 
mS 
,F 


1OO< : 

ria h, 

Tac TW z 
L bew 


biew ..1.9.4 


riew owl 6 Sae 


sonst. 


RaOps i 


wooest ,. Lida 


co HTsOn Ov agate TE CA 


ATAAHA Ui AavIa = ** ae 


Sag 26D 





Si.) good 


‘iS 


mp spe b BT stat 


\y 
rp 
ok 


S13 2] 7d 


iss teegonl 
lacataunad 


Letxseubai 


laixteubut - 


ipiasagont 


It Sagat - 


” 
a  nlnnennnnnnneinnemneitintnl 























a 
- ~ = 
- . a * “~*, 
“ > = = 
I wICS : 
o 

- + ae ~ oes 
eed 





garrol her ee apie to awot 
wort nodiend te aw. 7 

| owed 2: awa’ : 

sve ~ LAO Eokwegat ; 
wotnonbe, ‘20 wie ee 

bat ve baal’ entbennd 
2 aacubant pabbl bef , 


tj Sat Lio tstaeqaa’ 8 
: APA Ae ee soaxe? — : 








is 


TABLES Lies (cont'd) 





Character Type 
source of of 
Waste Treatment 
Town) Of Hort See reneren Domestic lagoon 
Sherritt Gordon Mines Ltd. Industrial lagoon 
Dow Chemicals iF alel Industrial lagoon, dis- 
posal well 
Redwater Imperial Fertilizer industrial lagoon 
Town of Redwater Domestic lagoon 
Redwater Imperial Oil Gas Plant Industrial lagoon 
Waskatenau Domestic lagoon 
Chemcell - Duvernay tecustpia. no, treatment 
Elk Point Domestic lagoon 
Canddtan salt Cos btd. Industrial no treatment 


ST - secondary treatment 
AL - aerated lagoons 


A.P.I. - API separator (American Petroleum Institute) 
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TABLE III 


EFFLUENT CHARACTERISTICS OF EDMONTON #5 SEWAGE TREATMENT PLANT 


1970 Flow Ss B.O.D. ee Phenols Coliforms ee 
mon. M.G.D. ibs/day 1bs/day Haya mg/l MPN/100 a 
neeeGe o6 “Tets03° 14" 256 6 18 8.1 = 
20) 830 -25peu2 1-16 2eee 6; 851. 15 24 ix 31.0 
Sees 99) 26,945 00 11,583 4 16 3.3 26.0 
4 34.36 29,066 14,470 5 11 ss, 34.2 
5 30.74 61,509 51,688 3 19 0.1 = 
6 32.06 63,895 72,606 = 96 7.9 40.7 
7 41.47 86,261 79,192 - = = - 

8 34.38 56,181 55,610 38 85 - 44.4 
guns 7099" “S4" 765°" 50,924 9 64 = 25.9 
10 32.38 48,209 43,981 6 24 = 25.1 
11 30.40 /23,528 23,364 9 15 = 24.8 
Pome Gens). | 21725 a G18 14 17 = 26.3 
EFFLUENT CHARACTERISTICS OF EDMONTON #3 SEWAGE TREATMENT PLANT 
i BE7L 3,518 4,040 22 42 2 = 
2 3.85 3,224 3,830 47 58 = 18.0 
3 4.47 6,469 3,623 26 43 = 19.0 
4 4.52 4,826 3,797 17 27 = 13.7 
5 4.72 4,006 4,617 20 30 - - 

6 4.36 2,954 3,880 30 25 = va Pes 
7 4.98 3,820 3,963 = 16 2 2 
g 4.59 3,044 2,906 18 30 = 22.4 
9 3.89 3,204 3,522 25 29 = 16.9 
10 3.50 2,829 3,638 14 24 = 20.9 
ll 3.08 2,629 3,696 26 26 = 2379 
1 3.25 2,514 4,046 37 34 = 2604. 


15 
is the major source affecting the river for a considerable 
distance downstream. Treated domestic waste is discharged 
into the river from the number 5 and number 3 (until 
December 1971) sewage treatment plants. Waste character- 
istics are ehimmageeee ite lables Li 

During the summer high flows the sewage from sewage 
lagoons located approximately 16 miles downstream from the 
Low bevel «Bridge is®dischargéd into the river. ~These 
lagoons store wastes from three packing houses, and house- 
holds@insthesnortheastern section of the cityeas well as 
from the Town of Sherwood Park. 

The Edmonton industrial complex has at times been 
neeroyei bie for odor problems found downstream from 
Edmonton, because of its phenol containing wastes. 

Farther downstream there are waste flows from the 
Town of Fort Saskatchewan (discharging domestic wastes) 
and SherrittGordon Mines Ltd. and Dow Chemicals Ltd., 
which discharge industrial wastes. The wastes from Dow 
Chemicals is of special importance to our study, because 
of its high chioride content, which could be traced a long 
way downstream and was used in this study to verify our 
calculations of variations in concentratim over the river 


cross-section. 
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ee eee 
THEORY OF DISPERSION IN NATURAL STREAMS 


PLD eLOduUCL ION. 

With the growth of industry and development of new 
technological processes the number of potentional pollu- 
tion sources increases as does the total concentration of 
pollutants in natural streams. There are limits, however, 
on how much waste streams can assimilate, and those 
limits are determined not only by the amount of flow, but 
also by the individual dispersion characteristic of every 
particular stream. The dispersion characteristics of 
saeetusveul streams vary greatly from stream to stream. 

A pollutant introduced into the river tends to 
spread out and become more diluted because of mixing pro- 
cesses fundamental to fluid mechanics and the dispersion 
characteristics of the river. 

In this study, an attempt has been made to study the 
mixing process with the use of dye and to predict cross- 
sectional variations of pollutants in the North Saskatche- 


wan River below the Town of Fort Saskatchewan. 
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The use of dyes in stream diffusion studies has been 
recognized for some time. Some detailed techniques concer- 
ning the use of dyes in measuring some of the physical pa- 
rameters, that are involved in the stream dispersion, have 
been developed. There are variations in the techniques pre- 
peneed by different authors, but all of the derivations 
have their fundamental assumption based on a finite element 
Seen using an elemental cube, as shown in Figure 3. 

The concentration of any dissolved substance in a gi- 
ven liquid in an equilibrium state is constant throughout 
chieswhole body of thesliquid. If the concentration is higher 
in one part of the liquid as compared with another,then, over 
a period of time, there is an increase in concentration in 
that part of the liquid that has a lower concentration and a 
simultaneous decrease in concentration in the part of the li- 
quid with the higher concentration. This phenomena can be ex- 
plained by an equal and random transport of liquid particles 
with their contained dissolved materials from one part of the 
body of the liquid to another. An assumption can then be made, 
that the rate of transfer of anymaterial in any direction is 
proportional to the variation of concentration in that direc- 
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Elemental Cube. 


FIGURE 3 : 





iS) 
justification of such an assumption can be given on the 
basis of the kinetic theory of gases and liquids for 
molecular diffusion (17). 
The infinitesimal cube as shown in Figure 3, has 
Svaes dx, dy, dz. tet “us consider, thats through the side 
AOEF, with the area dy dz, a quantity of substance equal 


cO 


a Oveczmat 


enters the cube during the interval of time dt. Through 


the side CBGD a quantity equal to 


ac 
DS Io xcax dy eka seks 
is leaving. The-minus sign is used to indicate that the 


movement occurs in the direction of decreasing concentra~ 
tion. The constant of proportionality, Dos is, Called the 
coefficient of diffusion for molecular diffusion and simi- 


larly is termed for turbulent dispersion phenomena. Because 
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20 
the difference between the quantity entering and leaving 
is equal to 

326 
D.— dx dy dz dt 
x 
OL 
x 
in which Dd. isechnesceoeceficient lof diffusion. Wifhertotal 
excess of material entering the cube equals the change in 
concentration within the cube per unit time 
2 ”, 
fap, 2f+p 284 0,28 (2) 
at emo se2 Y 9 2 ee 
a y Zz 
Ben etiixing sin Natural Stream. 

The mixing in natural streams is a very complicated 
proccess. s6l.C is known that 1t is a three-dimensional pro 
cess, in which vertical, lateral, and longitudinal mixing 
is determined by the geometry of the stream. Probably the 
most important part of this process is longitudinal mixing, 
which occurs as the consequence of the lateral velocity 
gradients. The surface velocity across the stream varies 
inedicect. relation ito the depth,-so that the deepest part 
of the channel has the greatest surface velocity and the 


SLeLuctiwObeLedecer dye: ane this» pPOrcion Of the echianned 


moves faster than in other parts of the cross-section. 
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In most natural streams the flow is turbulent and 
diffusion takes place in a vertical direction between sub- 
sequent longitudinal layers of water which have different 
velocities. Because of the relatively small depth of 
streams this Pacehoatl mixing! takes place much faster than 
longitudinal and transverse mixing. Vertical mixing can be 
taken into consideration only in cases in which the efflu- 
ent is soluble or forms a fine dispersion. Oily sub- 
stances with a specific weight less than that of water 
Wi PenOtemisin thesvertical direction, but will stay and 
Spread on the water surface. 

TiCmenICOmCOnSticuchtmOre the Mixing process , Lic 
ereveres Mixing, lsspLlimaralyedue to secondary currents 
existing in every stream. The three patterns of secondary 
currents are dependent on directional changes in the 
channel. Looking downstream, when the channel turns to: the 
right the mass of water tends to follow a counterclockwise 
spiral, which often persists behind the curve into the 
straight part of the channel. A similar pattern in the 
opposite direction is known to be formed in left turn 
curves of the channel. In both cases strong, spiral-like 
patterns of flow cause lateral mixing by bringing the water 


from one side to the other. 
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22 
In straight stretches of the river two spirals are 
usually formed with opposite directions of motion joining 
in the middle like a gear (15). 
The spiral-like flow in curves is much stronger than 
thatein) straight channels and therefore complete mixing 
could be expected to be achieved much sooner in curved 


channels. 


4. Basic Equation for Dispersion in Natural Stream. 

Using the same concept as applied for molecular 
diffusion, we can place the infinitesimal cube, as shown 
ine gure 93 sandwdesert-bcasineparte2. Of this»chapter,sinto 
the moving body of water, in our case into a natural stream, 
Similarly, mo can bemshowne that «the change tin concentration 


is equal to 
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This equation governs the dispersion of a slug of dye or 
any readily dispersable substance dumped into any turbulent 
stream. 

If there are no tributaries nor loss of water, the 


governing equation is more commonly used in the foilowing 












ce | = _ 7 ae 

oe alaciqn ov) weciy ont 25 Sqienan aietivednie ice 7 

galasat ‘om 2190220 Ge Peete pists hatte ig Saal 
Ai) tee0 i erat bhi ithe 


rewts 2 P. 4 ' . Walt wit ile-_ eto eee 


0 





' o colt Wis Bietaado WOR eee a eee 
apa WuSLioe SOO bet Seq Set biwap 7 
; . slananip 


a 
+, bn 
~ lie 


y ie Ci - * : Sd hot Cl Jt onl a) 


{Gus £6 cotvios Sube Sie eae - 


& = 
v| - i 
4a Ww mye 

} a BOAS aa 

* 

a » 

os ~ pA <a 7 = 
46. j _ ee | ©) el A} ss 


4a ah 20 ee a Io ore tena se a 


o <n waa ba? mete: 








_ 


ea Fa 


23 


form: 
BESET ni, fein OP RGASI hd cS aca [D =) es oe Sy 
at x dy 92 9x Ox ava Oy, 
3 ac 
+ [D,] (3) 
az az 
where 
c = instantaneous local concentration of the dye, 
in ppb by volume. 
Sonny pce stot lal Ccoordimatres in Longitudinal, 


lateral and vertical direction respectively. 
u, V, W = instantaneous velocity components in 
x, ¥, Z Girections respectively. 
Do, eee Dd, = turbulent dispersion coefficients in 
ft“/sec, inex ,ey, 2 directions respectively. 
t = time since the injection of the dye, in seconds. 
The first term on the left hand side of equation 
(3) is the local change in concentration with time; the 
next three terms are the change in concentration due to 
convection by the instantaneous velocity components. The 
right hand side represents change due to turbulent dis- 
persion. Molecular diffusion in natural streams is 


relatively small and is therefore neglected. 
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The solution of this equation is very difficult because 
Oreits Variable coerfLicients. In natural conditions varia- 
tions in u have the greatest effect on the amount of mixing. 
The parameters v and w are small or often zero. At the 
present time no gets on) SOlULTON Ofeccuatron (3) "exists. 
Various workers have devised some approximate solutions, 
sical based on a one-dimensional dispersion model, by con- 
Sidering variations in the mean flow velocity in the cross- 
section as the primary cause of dispersion. 

The following equation was originally developed by 


Taylor (5) for this one-dimensional model: 





where 
lie=eaverage velocity of flow 
DS = constant longitudinal dispersion coefficient 
(D, = D,). 
Several assumptions have been made in deriving this 
equation: 


Ie tWatereislincompressible! and! itsevélocity ulis 


time independent. 
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25 
gee needensaty of the.solute issequal.to. that: of 
water. 
Jame Leche Olcbersi2Onncocificient Dy, is constant. 
4. The mixing is uniform and the mean concentration 


Cc depends ColysOnexeand tt. 


thessolution BotLeecquation. (4).is then: 


(x - ut) *, 
4D, t 


M 
Circa) ei STEELE exp [- G5) 
A(4nD,t)” 5 


whererc(x 70) == 075 for x>0. 


The boundary conditions which must be satisfied are: 


pe A) clint medx —eM forall vyatues of £ 
c(m,t) = 0 t>0 
che >0 as tro 
ox 


where M = total volume of dye injected, 
A = effective cross-sectional area of flow. 
Equation (5) gives only a theoretical distribution 
Of pollutant or dye concentrations ;<Ets applicability to 
natural streams has often been questioned. 
In the practical measurements of dye dispersion, an 


unknown amount of dye is lost due to photochemical decay 
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26 
and benthic adsorption by vegetation and suspended solids 
present in the stream. Variable amounts of dye may be 
temporarily detained in and later released from dead zones, 
which are usually present both in lateral and vertical 
directions along a ane orga! stream. 

Eouation oymaoesenot, take account of these two impor t— 

ant factors and therefore its applicability for practical 
purposes remains restricted. Factors K and K in the 


0 


following equation better express the above described 





effects. 
1 M eee x 
CS = ee ae ——] (1-K)—) (6) 
K §A(47D_t) a Abe Ve 
L L 
where 


K = regional dispersion factor 


Ky = coefficient for the loss of the dye in the 
reach 
V = the average velocity of flow 


The value of the regional dispersion factor K is 
expected to vary from stream to stream. Values Gt ok muss 
be found by the trial andj error process of comparing the 


computed values with those observed. 
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27 
The value of Ko is expected to depend upon the 
following parameters: 
Dy = the longitudinal dispersion coefficient, 
Peeacicmcopswiatheor tlow, 
S =) the oe to peak arrival of dye concentration, 
x = the reach length, 
Pemeticmenasacterictic. lengrh to which the Pateral 
eeoevertical mixing, 1s complete, 
Ki peciieomenerricient tor the loss of dye due to the 
decay and absorption. 


A dimensional analysis shows that: 


B H L iL 
Ko =f ( a ' ae, ~F =, cache K,) (7) 
Dt Whe, 8, L VB 








Ticwcenpara cal Lormudast or) thesdistancesin is: 





16 D 
L = = (8) 
0.05V 
and for coefficient Ky 
H 
Log (2) 
Ky = =f ME + ——2- + I (9) 
ite 20 pe Logtby—2) 
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In the equation (9) 


Peet actemeaccouUncIngqe1or the detention of dye in lateral 


dead zones, 


log (22) 


ct 


E 


Nh 


Comcmi acco eacCOunmine fOr the detention Of dye in verti- 
cal dead zones, and Ky tena cOertici1ent detained fora given 
estate (ee: the loss of dye due to decay and absorption (in 
our study Ky =e) Peo eSoeCOC LE TCU Chom lS stake as, ca =COn= 


Stant, but it is rather time and concentration dependent. 


The resulting final solution of equation (4) is 


then: 
J M feeVE) Cc. 
c(x,t) == Tee SS l= J 
K A(4nD_t) MD 
L L 
x B xX log (2 Ly x 
Le eee Cae a Ky ave (10%) 
16D,t 4vt log (D, —B) 
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29 
This solution of the dispersion equation satisfies all 
the wnecessary initial and boundary conditions required in 
EnewCa5ceOLrea Naturalsstream. It may be vused to calculate 
the dispersion coefficient for different reaches of the 
stream where time Benet Verian curves are available, and 
LOY a given dispersion coefficient it can be used to deter- 


mine time-concentration curves. 


J eelev lous. oLudles, 

The first study of dispersion phenomena was carried 
Outre bysTaylor (5) eusing the concept of longitudinal dis- 
persion as a one-dimensional process. He limited his 
Scud vecCOmanL LOW au tOongestrailgic pipe Lor which be ob- 


tained dispersion coefficients. 


Oo 
1 


MOAR ey 1oRY (zis) 


where, 

D —- dispersion coefficient 

a —- the pipe radius 

U* - the shear velocity 

Applying the same concept as Taylor, Elder developed 
a formula for an infinitely wide two-dimensional channel and 


obtained 
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IY Ss Bagels) Goh hube 2) 


where d = depth of flow. 

Tay LOS concept, derived only for the flow in 
pipes, Nas become widely accepted and it took further 
studies to prove that dispersion coefficients for natural 
streams are in fact of much higher values. Engineering 
practice since then has been to assert that a coefficient 
exists, but because of specific stream characteristics it 
must be obtained only by experiment. 

Various investigators have derived empirical dif- 
fusion formulas applying to the flow in pipes, in open 


channels under laboratory conditions, and in natural 


StLeans sew OLPougestudy, only Lomuulas for natural streams 
aveeoLeconcern (9). -errhey are: 
Glover Ds, =95007U* 9H . (13) 
preereconntn. 12 ylasdue pele (14) 
L u* 
a i 1 f 
Fischer D = -= 4q'(z)dz 4 ——-— __ dzg q' (z)dz 
A BE d(z) (15) 
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= overall bed shear velocity 

hydraulic depth of flow 

mean velocity of flow in direction of flow 
area of flow cross-section 
= depth integrated velocity 
= transverse turbulent mixing coefficient 
transverse cartesian coordinate 

Wlath of Channel 


depth of flow 


Dimensional Mathematical Model of Dispersion. 
one-dimensional model is expected to be valid 
ross-sections at a considerable distance down- 
m the point of injection, where the mixing is 
oubeminvtornmee Nearmthe point Of inyection the 
incomplete especially in y and z directions. 
womterms On @heyright (‘handisadeyor eq.) (3) can= 
lected, but the convective terms due to mean 


omponents in y and z direction can, as they are 


usually small compared with u term. The equation can 


thus be re 


written as follows: 








Ther t~nitial conditions are: 
CH yy Zee FOREX 0 py 0 pear. 


The boundary conditions are: 
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Glory, 2, cj = 20 for all values of °t 
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Assuming D,. = Dy 


velocityeyv. for the seach, the solution “of eq.) (16) is 


and the mean velocity u = 
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(20) 
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This equation satisfies the initial and boundary 
fondue Onsmctvenes) @eo suai )yuat)),, (20)) = butett. doOeaenot 
satisfy the boundary condition given by eq. (18). However, 


eq. (21) does satisfy the similar boundary condition: 


co foe) foe) 


J J eec(x 7 y2,b)1dx,ay dz.= M (22) 
0 
LGrealie values Of € 


Incorporating the preceeding, eq. (18) can be hypothetically 


approximated as: 
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eo By eee 
= wh i Cera, Zit) ax dy az 
foe} co H 
<2 RV omy» RVSatnegi sks ctr ae 
0 B/2 
ap -B/2 a e(opy zt) dxedy edz (23) 
Eq. (20) indicates that for natural flow conditions 


in a stream with finite values of B and H, the amount of 
dye recovered would be less than the total amount of dye 
injected. The amount of dye lost or detained in dead 
zones laterally and vertically is assumed to be given by 
BhewcilLcerquantreies fone the right hand side .of eq. (21). 
Eq. (21) still does not take into the account the 
amount of dye which may be lost due to adsorption or 
decay, and therefore could only be valid for a short time 
after injection of dye, when dye losses are relatively 
small. For distances further downstream some modification 
of eq. (21) may be necessary to account for such losses. 
Eq. (21), known as the three-dimensional mathemati- 
cal model of dispersion, has been used in this study to 
determine dispersion of dye in the reach of river from Dow 


Chemical's outfall at Fort Saskatchewan to the government 
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ae 
sampling station at Vinca, where measured time concentra- 
tion curves were available. 
To use this equation for three-dimensional mathemati- 
cal determination of dispersion it is necessary to know the 


coefficients Dor D D_. The value of Dy. is assumed to be 


Yaa Z 


equa tomD s-sthe Longitudinal dispersion coefficient 


L 
determined from one-dimensional model of dispersion. How- 


ever, dD. may not necessarily be equal to D Oh serene 


L’ 
Taeetiay a UessOrectener sxsor t, it would approximate 1. 
The lateral dispersion coefficient, Bho is a £unc— 
Elon ob the vertical, dispersion coefficient, D, the top 
Width, B; the depth of flow, H; the average velocity of 


flow in the reach, V; and the mean velocity of flow at the 


sampling station, Ve dimensional analysis shows that: 


= ve B 
Dee eras >.) _ (24) 


where Ape Bs 2; V i and Q = discharge, in cubic feet per 


Pp 
sec. 
The following is the empirical equation for DY used in this 


study: 


Dees (7) D (25) 


_ ee 
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36 
The value of DF which gives the best fit between the 
computed and the measured time concentration curves is 
accepted as an appropriate value of dD, atsrhat’ Location of 
the stream. For practical calculations the value of Dd, is 
usually assumed. In our case the value of dD, was assumed 


5 


to be 1.52 x 10. Pee ced. 
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CHAPTER: IV 
DISPERSION COEFFICIENT DETERMINATION 


1. Mathematical Analysis. 

The eq. (10), as Statedaine Chaprer 1 wseithessolu— 
tion of the basic equation for the one-dimensional mathe- 
matical model of dispersion.gs ite 1sanecessary (topmeasure 
several physical parameters of the particular stream to 
complete the solution. The most important factor is the 


Jongitudinalsdaspersionecosctiacient) D As explained in 


L’ 
Chapltermm EE wpanteoue toiesalspersionscoctficient must be 
obtained by experiment. There are several methods for 
determining thes dispersions, coefficient, but in this case 
eq. (15) and calculations based on physical data of the 
river seemed to be most convenient. 

TOmUSemed MELO meLObethe Gisperston) Cocrtrticient cal= 


culation, the integrals must be replaced by the summation, 


Lesulting eineeiew:Ollowing ~expresston: 


Cz q'4,)] (26) 


in which q' = > (d, + 4, 
i 
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u, = mean velocity in i-th vertical slice u' = 


Ufa u , where ; 

Uf =amnean velocity of flow within entire cross- 
section, 

a. = depth at the begunnang) ofthe wath vertical 
slice, 

a On 23 d; U* = transfer coefficient between i-th 

angie) ethivertical slice, 

Ne Number iotavestucalsrsilaices (notyless»than 20), 

A = area of flow cross-section, 

U* = the overall shear velocity. 


The shear velocity may be calculated 


U* = (gr Shes where 

g = the acceleration of gravity 

©” = the’ hydraulic radius, 

S = the slope of energy gradient. 
2. Field Measurements at Vinca. 


The equation (26) indicates the necessity of deter- 
mining the velocity of flow and its variation across the 
river as well as determination of other physical data. For 


this reason, the cross-section was measured and divided 
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39 
into sections 20 feet wide. Holes were drilled through the 
ice every 20 feet to measure depth, ice thickness and veloc- 
city of flow. The ice thickness did not vary significantly, 
being between 15 and 17 inches. The number of stages for 
flow velocity meseucemente in each hole depended on the 
depth. The river is quite shallow and so three stages were 
taken in only three or four holes. An hydraulic current 
meter with several exchangeable propellers was used for the 
actual flow velocity measurements. The choice of various 
propellers suitable for various ranges of velocities en- 
abled greater accuracy achieved. 

The data was then summarized and used to calculate 
the fiesereion coefficient De (Table B of Appendix). The 
correctness of the field measurements of flow was checked 
with the government record from the gauging station at the 
Low Level Bridge in Edmonton. The calculated flow of 
2089 c.f.s. is very close to the government recorded flow 
OLe2zi205Cl£.S. (Since there are no significant tributaries 
between Edmonton and Vinca Bridge, and the small tributar- 
les are frozen during the winter months, the difference is 
close to 1.5 percent. 

The characteristic cross-section, with velocity con- 


tours interpolated from the velocities obtained, is shown 
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40 
in Figure 4. In this figure the sampling points used for 
the dye dispersion study at a later date are also shown. 

The dispersion coefficient calculation is based on 
Pisher!secomulan(!0) equation (15). The resulting value 
of the dispersion Peri cient aS 9 / OA 0% Om sec. Despite 
great care during the measuring procedure, this coeffici- 
aoe PSenOuw Stiltable for sehe time-concentration curves cal— 
culation for reasons described later in Chapter VI. The 
dispersion coefficient based on Glover's formula (9) equa- 
GLOuss La ewe ene value 230 fel sec: appeared to be more 


Suitable for the river stretch being investigated. 
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TABLE LV 


RHODAMINE W. T. CHARACTERISTICS 


Dye Color 

Chemical Formula 

Form Used 

Specific Gravity 

Minimum Detectable 
Concentration 

Photochemical Decay 

Physical Adsorption 

DeLiusiv icy 


Biodegradability 


Basic Violet 10 

Co gH3,N503C1 

20% solution (by weight 

Ug de 

0.000046 mg/l in distilled water 
0.000066 mg/l with background present 
50% concentration reduction 

moderate 


fair 


ies be 
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CHAPTER V 
DY EV RACER STUDY 


ie cevyLlous Studies. 

Several studies have been carried out on the North 
Saskatchewan River, mainly to determine flow velocity in 
ay. PeacieseOne tne river, (Re. D. Cameron (1) and 
DRBCOEGitalre (29) == [he results” of mHOee studies do not 
entirely reflect the present situation, because the ice 
formation varies from year to year. However, for rough 
estimation of the time of travel of dye, they provide 


valuable data. 


2. Measuring Equipment and Fluorescent Dye Used. 

In this study the dye was injected instanteously 
as a Slug. Based on experience reported in previous studies 
where either 25% fluorescein or 40% rhodamine dye were used, 
Rhodamine W. T. dye was chosen as the most eee 
Table IV lists some important data concerning this dye. 

A G. K. Turner Associates Model III fluorometer was 
used for actual concentration measurements. The fluoro- 
meter 1S basically an optical bridge, analogous to a 


Wheatstone bridge used in measuring electrical resistance. 
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The optical bridge measures the difference between the 
light emitted by the sample, which is proportional to the 
concentration of fluorescent material in the sample, and 
that from a separate calibrated light path. The quantity 
of light required La mMiceecali brated lightwpath to balance 
with that from sample is automatically indicated on a 
dial marked in 100 equal divisions. The fiuorometer is 
equipped with four openings or appertures which emit eit 
to the sample in the approximate ratio 1:3 : 10:30. ‘This 
provides a means for examining a wide range of concentra- 
tions without changing the filter arrangements. The pur- 
Bose OmecOlour f£i2tters: in a fluorometer is to limit the 
i aa the photomultiplier, as far as possible, to 
that emitted by the dye. 

ie actual concentration of dyé in a sample is de- 
termined by multiplying the dial readings by a calibration 
constant. sihe calibration constant 1S a ratio of the known 
concentration of the dye C to the reading obtained for this 
concentration D. Care must be exercised that the sample 
reading as well as the calibration constant are obtained 
within the same range. 

All measurements were made using Lhesgencrad purpose 


ultraviolet lamp and the standard cuvette holder door 
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supplied with the instrument. Cuvettes were borosilicate 


glass 12mm x 75 mm tubes. The zeroing cuvette was a round 
black 12 mm x 75 mm tube simulating a non-fluorescent 


sample. 


3. Dye Dispersion Measurements. 

As mentioned earlier, the suspected incomplete mixing 
of wastes from Fort Saskatchewan to Vinca established the 
need for this study in that’ locality. The main source of 
detectable chlorides, Dow Chemical Ltd. plant, determined 
also the first injection point for the dye. The sampling 


point for the passing dye cloud was then established at 
the Government sampling station approximately 2 miles up- 


stream from Vinca Bridge. The net travel distance as 
measured along the river centreline is 70,600 feet. 

In the second run the degree of mixing between 
Edmonton and Fort Saskatchewan was examined. The release 
of dye was below the Edmonton water treatment plant, 
sampling was at the Fort Saskatchewan Bridge. 

Listhicsei CString 0709 mol 2045 Rnodamine Wea te 
was released into a manhole located about 100 feet from 
the river on a sewer outlet line from Dow Chemical Ltd. 

It was estimated that the dye would reach the river in one 


or two minutes because of a relatively high flow in the 
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46 
sewer. The dye was released at 0:00 February 4th and 
measured in the cross-section 13 1/2 miles downstream 
between 11:30 A.M. and 3:40 P.M. on the same day. Time 
of travel as estimated in the previous studies was within 
the time range established by peak concentrations which 
arrived between 12:20 P.M. and 2:45 P.M. and were measured 
at six points across the river cross-section. 

The actual sampling was carried out from six holes 
drilled through the ice. To ensure that the passing cloud 
of dye had arrived, a fluorometer was set for continuous 
operation on the hole where the strongest flow in the 
river was detected in previous measurements. Low aixz temp- 
eratures caused freezing of the pump supplying the eieees 
meter, but fortunately only after it had shown a sharp in- 
crease in concentration. The passage was then established 
from bottle samples taken at intervals of 10 to 15 minutes. 
The resulting curves are shown in Figures Al to A6 of 
Appendix. 

The second dye run was carried out between Edmonton 
and Fort Saskatchewan on Thursday, February 17th. 34.67 lbs. 
of 20% soiution of Rhodamine W. T. was released at 10:20 A.M. 
just below the intake of Edmonton Water Treatment Dillant. 


The actual sampling took place on the following day between 
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47 
10:30A.M. and 3:40 P.M. below Fort Saskatchewan Bridge 
from 4 holes in the ice. The time-concentration curve 
obtained (Figure 5) showed that the mixing from Edmonton 
to Fort Saskatchewan was sufficient to obtain a relatively 


uniform concentration across the river cross-section. 
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CHAPTER VI 


APPLICATION OF ONE- AND THREE-DIMENSIONAL 


MODELS OF DISPERSION 


1. One-Dimensional Model of Dispersion. 

Equation (10) is the solution of the one-dimensional 
dispersion equation that satisfies all the NecesscacyveLnati— 
al and boundary conditions required in the case of a natural 
stream. The concentrations for time-concentration curves 
were obtained from the first dye run, described in Chapter 
V, part 3. Using the dispersion coefficients determined 
in Chapter IV, an attempt has been made to calculate time 
concentration curves by equation (10). Since the one-dim- 
ensional model of dispersion applies only to mean concen- 
trations in the cross-section, the mean concentrations had 
to be calculated for each time interval from our measured 
data and existing flow. This "dye" curve was then compared 
with the curves produced by the equation (10) using both 
dispersion coefficients. The data for all curves are shown 
in Table V and plotted in Figure 6. 

it is apparent from Figure 6, that the calculated 
curve uSing the dispersion coefficient 7040.5 ft.*/sec. 


(Chapter IV) is too flat as compared to the "dye" curve. 
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50 
It suggests two important facts: 1) the dispersion coef- 
ficient used is too large, and 2) the one-dimensional dis- 
persion model does not fit very well the situation where 
mixing within the cross-section is not uniform. 

For the cee possibility a coefficient based on 
Glover's formula (9) equation (13) has been used and the 
aiseanee obtained (shown also in the Table V and Figure 6) 
Show closer resemblance with the "dye" curve. 

Despite the fact that the expression used takes into 
account the detention of dye in dead zones, the tailing on 
the measured curve is much longer. Using the trial and 
error methods dispersion coefficient has been assumed 
DS = 980 fe, Veec, The resulting curve, Shown in Figure 6, 
gives the best approximation to the measured data. 

This one-dimensional model, even if it produces a 
close approximation of the measured curve, is valid only 
for a cross-section a considerable distance downstream from 


the point of injection, where mixing is more uniform. 


2. Three-Dimensional Model of Dispersion. 
Equation (21) represents the solution of the three- 
dimensional dispersion equation. The same expression for 


the loss of dye as in the one-dimensional model has been 
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TABLE V 


TIME-CONCENTRATION CURVES DATA INDICATING THE PASSAGE OF 
DYE FAT, VINCA 


(BASED ON ONE-DIMENSIONAL MODEL OF DISPERSION) 
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Sine —— riebramenntnneete tiie Steet eel 
D.=7040.5 £t“/s D.=230 £t*/s D.=980 £t7/s 
L L L 
hours ppb ppb ppb ppb 
ete) OR ZO Om7-00 On027,9 05320 
11:40 (0) areah 
jE 210) OMS 49 Ona 0 OO 7eks Bers 
e200 0.404 
OE0 ORD 5.6 ON oi Oma 0.484 
DREAD On G26 
0:30 Om o73 05745 Ons23 0.566 
0:40 OOo 
Cs53 0.760 Oa 60 0.548 0.644 
205 Oxy Dy, O03 7.60 0.623 ORG TE 
Degas) 0.716 
13.0 aes 32 0.741 Ons 2 Oa a2 
1:40 One 21 
AG SKY) Onu/L6 O34 Oe70>5 OR 09 
ZG She JASN 
2rl5 ae ley, ONL Od 2 0.681 
BEd Om6o7 
2240 0.624 OnG93 OR esa. 0.624 







: 


as 

- Ati | : . _ 7 

a3 DATTADIOMY (Ages Savi? iC aite 
a a v te 


1 


V Stra? 


i 


ors Pa vA ied 4) 


0 * i 
io sse ae 1 LITO CAND Seat ae oxvesg) 





Sy, 








a aoe ee 00:2 Of :T 00:1 0€ 10 00:21 Of :T1 
Bf cree nara (atari | ieee |e meee ieee | Fellas al |i 2 aren ete el ie, ol ie eee IS ae! et). ace 
0 
a 
ee 
ee 0° 086 =1q POIeINOT eS) ae ne A Lee 
8/233 0°0€Z = per,e{nope9 —90——-¢ —— ye 
8/73 GC OV0Le= (aepea et home) ota rs 
peianseapy vi : 
/ ce 
/ L 
i a 
x ys oy 


—— 
a 
~~ 
aX 
~~ 


~~ S x 
eth Oe _ a fo __e a —vV 
anes ——o¢ — as 4 Sha 


Gol oar avs Pee ree) dT akq JO aSesseg SUTQYOTpUT SeAIND UoTIeAQUaDUOD-oeWT] : 9 aaNdId 


qdd ut wvotjeirquaou0g 


. [o at wet | nor =f OsotL 
= 
-* 
(nawetocacen! Jy » Gateelat) ani hp a © egetsa® go =3e2i10h ? Shed on Jee See 


ge | 





53 
incorporated. It takes into account the loss of dye due 
to decay, adsorption and loss in lateral and vertical dead 
zones. 

Nearmacnehpoinc or injection tne mixing of dye in y 
and=2Z, on lateral and vertical directions, is not complete. 
This isjalso true in our case, even though our profile 
is located 13-1/2 miles downstream from the point of in- 
jection. The observed character of flow in this cross- 
section together with the measured time-concentration 
curves suggests that the centre of the axial system is 
200 feet from the south bank. All the samples used for 
the time-concentration curve were taken from below the ce; 
close to the = tce-water boundary, so the yz plane has been 
set also on this water level. The z dimension has there- 
fore been considered to be 0 in the following calculations. 
EhiemvalLuce of W762 08 ap has been assigned to the disper- 
Sion coefficient Dy in z direction as described in Chapter 
Iii, part 6. The time-concentration curves have then been 
calculated at the same locations as the measured ones. The 
Spatial simulation of results obtained as shown in Figure 7 
as compared with measured ones, Figure 8. 

For better comparison the calculated vs measured 


time-concentration curves at appropriate locations are 
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pewetnoqtes = 9 


peazanseow = W 
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by 
shown in Appendix A. Table VI shows the calculated results. 
Again, as in case of the one-dimensional model, there 
is considerably larger tailing on the measured curve as com- 
pared with the calculated one. The peak values on both 
curves are approximately the same. For this case of a 
three-dimensional model of dispersion the longitudinal 


dispersion coefficient based on Glover's formula, equation 


(13), has been used. 
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CHAPTER V Leb 
LABORATORY TESTS 


ipameeetroduction™. 

During the winter period 1971-1972. an, extensive 
Sampling program was ‘carried out both at Vinca and at Fort 
Saskatchewan inorder sto,obtain, anformation about. waste 
levels in the river water. Series of samples were taken 
across the river to determine the variation in concentra- 
tion. 

Various tests were run on the samples obtained, but 
only a few indicated the suspected incomplete mixing at 
the Vinca Bridge sampling station. Results from samples 
taken at Fort Saskatchewan Bridge confirmed fairly good 
mixing of wastes from Edmonton. There was no significant 
Variation in concentrations across the river at Fort Sask- 


atchewan Bridge. 


2. Laboratory Tests. 

The study area is located in the vicinity of Edmonton 
and the effects of its industrial and domestic wastes on 
the river water quality are felt far downstream. Several 
tests were chosen to indicate the level of various pollut- 


ants from Edmonton and from some sources in Fort Saskatchewan. 
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The best results were obtained by determining chloride 
concentrations. Chlorides occur in natural waters in 
widely varying concentrations. In upland waters chloride 
Contenteus usually quite low, riversand ground waters, 
however, show sometimes very high concentrations. [In our 
case, since the North Saskatchewan River rises in the 
easter Mountains as a mountain stream as do its tributaries, 
i omlasmamVcnvuelowicencentratlonvot chlorides S.) Asi tl des— 
cends into the plains and absorbs wastes from industries 
ancenunve tpalwiines sris=chloride content “rises, but ethe 
overall content of chlorides stays far below 250 mg/1, 
which is the recommended maximum level in a public water 
supply. 

It has been estimated that average daily excretion 
of chlorides is about 6 g per person or about 2,760 
kg/day from the City of Edmonton with its population of 
460,000, which would cause an increase in the mean con- 
CevcraG@ionwot 0 s5amg/1s, Glindustriest*also comteribute sto the 
Potro lmchlorwdesmcontent sy) somthateanotheresiqnificant i= 
crease is found below outlet from Dow Chemicals’ plant in 
Fort Saskatchewan. The Town of Fort Saskatchewan and 
Shermciute Gordon Mines Sip lant also tcontributermcomtme overall 


Contentmof chlorides tinethe waiver nn ourestudy tarea; but 
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60 
their contribution is small in comparison with that of Dow 
Chemicals Ltd. 

Another test applied in our study was for phosphates. 
Phosphates are added to the river via industrial waste 
discharges from alante using phosphates to control scaling 
in boilers, or via municipal waste discharges, where they 
ocette as result of human metabolic breakdown of proteins 
or from household detergents. They occur as orthophos- 
phates and polyphosphates, the latter being slowly trans- 
formed to orthophosphates. The speed of this transforma- 
tion is greatly influenced by temperature. 

An indicator related to the concentration of organic 
merece e the dissolved oxygen level, therefore the test 
for dissolved oxygen was chosen as another means to deter- 
mine distribution in the river cross-sections. 

Bacterial fauna in the process of decomposition uses 
oxygen from the water and the oxygen deficiency at any 
point is then a measure of the amount of organic matter 
previously decomposed in the stream. The ice cover, which 
prevents reoxygenation, favors this method. 

A similar test, but one that gives more accurate 
results, is the test for total carbon (TC) and total organ- 


ic carbon (TOC). This test was run by Federal Department 


O£ Health laboratory. 
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CHAPTER. VILLI 
CONCLUSIONS AND RECOMMENDATIONS 


l. General Conclusion. 

The goal of une study, as formulated in the first 
chapter, has been achieved. There are, however, several 
points which need to be refined and further verified, but 
in the short time period of ice cover on the river, fur- 
ther reduced by bad weather conditions, it was not possible 
to accomplish them by existing means. It is assumed that 
fUeerer studies will) be carried out on the river, either 
by the University of Alberta, Department of Civil Engineer- 
ing; ai by the Government of Alberta, Department of Envir- 
onment and that this study will be of some help in this con- 
tinuing program. iso the applicability, which is presently 
limited to low winter flow and perhaps only to this year's 


ice formation, will be extended. 


2. Determination of the Dispersion Coefficient. 


There are several methods used in this determina- 
tion, but the one used in this study was expected to give 
the most precise solution. Even though real care was 


excercised during the field measurements, the dispersion 
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62 
coefficient obtained was not suitable for our situation. 
The coefficient based on Glover's formula, however, seems 
to give more reasonable results, but the best approxima- 
tion of the measured data provides the coefficient 
De = 980 ree which was determined by trial and error 
method. Should there be any future study an attempt should 


be made to verify this coefficient for different flow con- 


ditions and possibly for various stretches of the river. 


3. Concentration Distribution in ChemRivereCross-section 
at Vinca. 

As is shown in Chapter VII, the best results, as 
far as the concentration distribution over the river cross- 
section at Vinca is concerned, are given bytthe test <for 
chlorides. These results, as shown in Appendix (Cy indi 
cate the highest concentration levershontthe vightespac "oF 
the river. The data measured and calculated for the dye 
tracer run indicate that the ratio of peak concentrations 
at various sampling points within the cross-section to the 
amount measured at the approximate position of Government 
sampling point is as follows: 

Ratios Mn Se dye released at Fort Saskatchewan 
(right bank) and measured at Vinca are given in the following 


table. 
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Distance from the right bank in feet: 


100 200 300 400 500 550 600 Mean 
Measured Cp) mee Cee Oe) 204) ae 00 0544 9 163 


Poet Ome 4 Ogee mee oO 6 3.04 Nib bee le 00", 0.61 | 1e84 


oA = concentration of the dye at the point n 
C = concentration at the Government sampling point 
nh = distance from the right bank in feet 


Looking at the plotted measured and calculated peak 
concentrations, Figure 9, we can see a very close resem- 
blance, which is in contrast with the above tabulated data. 
The difference between corresponding measured and calcu- 
lated values is very large because those ratios are based 
on concentrations (measured and calculated) at the point 
of Government sampling, where, unfortunately the differ- 
ence between measured and calculated values is largest. 

The ratio of mean concentration to the base value is of 
special importance, since it can be used to approximate the 
Mean CONCentration an Che profile on the basis of thesusual 
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4, Relation Between Values Obtained at the Government of 
Alberta, Department of Environment Sampling Point and Mean 
Concentration in the River. 

For dye applications the ratio of values measured at 
the government sampling point to various concentrations in 
the cross-section, as well as mean concentration, is given 
in the table in the previous paragraph of this chapter. 

One must itake Gntosaccountimany other factors’ to’ make 
practical use of those values. The approximate application 
to pollutants is possible only for chloride concentrations, 
and only after the mean concentration in the river above 
Fort Saskatchewan has been determined. Another important 
fact influencing the accuracy of this method is the dis- 
charge from Fort Saskatchewan lagoons and Sherritt Gordon 
plant. 

The use for other pollutants remains virtually 
impossible mainly for two reasons: 1) The main source of 
incompletely mixed chlorides is Dow Chemicals, whereas 
sources of other pollutants might be as far upstream as 
Edmonton or further and their mixing is much more complete. 
2) The level of their concentrations is too small for 


Practical work. 
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Se actical Application Of Dye Study Data. 

The two dye tracer runs as described earlier pro- 

ide valuable information about the degree of mixing in two 

river cross-sections. The three-dimensional model of dis- 
persion used to calculate concentrations throughout the 
cross-section could be also used to calculate amount of 
dye or pollutants released at Dow Chemicals. An attempt 
has been made to verify this calculations, but it appeared 
to be extremely difficult to estimate how much of the 
total amount of chlorides present is being brought by the 
river upstream. Because of continuous discharges of 
wastes a much simpler mixing equation can be used to serve 
the same purpose. Assuming that the concentration of 


pollutant remains constant throughout the day: 


QceH= G1oy + Ugly + 4303 + oe eee eet Gc, C2) 
where 

Q = river discharge 

© = mean concentration in the river 


Chit discharge of tributaries 


Gy = Concentration Of tributaries 
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67 
Available data from our own and government tests 
were used to show an example calculation of concentration 
for the discharge from Dow Chemicals. Unfortunately avail- 
able data did not quite correspond and effluent concentra- 
tion from Fort Saskatchewan lagoons was not available and 
hadietombe: Sstimated.as the! calculated value of 5550 mg/l 
eairile be considered to be in the same range as the measured 
concentration of 6250 mg/l obtained by government two days 


earlier. 


6. Recommendations. 

Bad weather and short time period did not allow more 
data to be obtained, which would provide a firmer basis for 
the conclusion. This fact itself suggests further verifica- 
tion of dispersion equation application and more detailed 
determination of the dispersion coefficient should be done 
for various discharges and flow conditions on the river. 
Further study of discharges and their concentrations is 
necessary in order to have a better picture of the overall 
pollution situation, which is now largely based on one or 
two grab samples analysis per month. 

The government sampling station, which takes samples 


about 150 feet from the left bank at Vinca, should make use 
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of the above mentioned coefficients in order to obtain mean 
concentrations in the river. The flow and dispersion pat- 
terms should be examined from time to time in order to up- 
date the information concerning the concentration distribu- 
fron, witcouwrs ea only for winter flows and most probably 


this year's ice formation. 
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CHEMICAL LABORATORY METHODS APPLIED 


Weelest tor Chlorides. 

Chlorides in water may be detected by means of 
volumetric procedures. In this case the Mohr method 
employing silver nitrate as a titrant and potassium chro- 
mate as indicator was used. 

As recommended in "Standard Methods" (18) 

0.0141 N solution of silver nitrate was used. This normal- 
ity of silver nitrate solution was chosen so that each 
milliliter is equivalent to 0.5 mg. of chloride ion. In 
the titration the chloride ion is precipitated as white 


Silver chloride 
oe =e 
Aguutacig = AgCl 


To detect the endpoint an indicator must be 
present to demonstrate the excess of Ag™ LONs ot Our Case 
potassium chromate was used as an AndacatorpepAsumchne con= 
centration of chloride ion approaches extinction, the silver 
ion together with chromate ion begins to form a reddish 


brown precipitate of silver chromate. 


2ag’ + CrO, = Ag,Cro 
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Presence of this precipitate is evidence that 
all the chloride has been precipitated. Since an excess 
Goas silver ion is needed to produce a visible amount OF 
silver chromate, the indicator error or blank must be 
determined and subtracted from all titrations. 

If the silver nitrate solution used for CitGa— 
tion is exactly 0.0141 N, the Galculation ron chlonides 
is»then: 

(m1 AgNO, = plank) 1.020.000 


Cle insti.) ei ae ania 6 
ml sample 


Pe Test for Phosphates. 
Phosphorus compounds occur in several varieties, 


the most common are: 


orthophosphates - trisodium phosphate Na,P0, 
- disodium phosphate Na HPO, 
- monosodium phosphate NaH.PO, 
- diammonium phosphate (NH, ) HPO, 


polyphosphates - sodium hexametaphosphate Na, (P03) ¢ 


- sodium tripolyphosphate NacP 3016 


- tetrasodium pyrophosphate Na ,P504 
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All the polyphosphates gradually hydrolyze in aqueous 
solution and revert to the ortho form. The rate of reversion 
is a function of temperature and increases rapldlyaaseche 
temperature approaches the powldinga pol nteweTheSratesis 
also increased by lowering the pH. Advantage is taken of 
both facts in determination of total phosphate concentra- 
tion. In our study the same method was used to determine 
orthophosphates and total phosphates. Polyphosphates were 
calculated as a difference between those two concentrations. 

The same colorimetric method was employed for both 
and the only difference was the pretreatment of samples for 
total phosphates by boiling for ninety minutes. By boiling 
Lie eee the polyphosphates are converted to orthophos- 
phates and as such they are measured. 

The stannous chloride method as described in "Stand= 
ard Methods' is based upon the fact that the phosphate ion~com- 
bines with ammonium molybdate under acid conditions and 


forms a complex compound known as ammonium phosphomolybdate. 


a + : i 
PO te 12 (NH,) MoO, + 24H = (NH, ) 3P0, 12Mo0. as 21NH, 
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The molybdenum contained in ammonium phosphomolybdate 
is readily reduced to produce a blue colored sol that is 
proportional to the amount of phosphate present. Excess 
ammonium molybdate is not reduced and therefore does not 
interfere. Asa reducing agent stannous chloride solution 


in glycerol was used because of certain advantages it has 


over other reducing agents. 


A+ 


(NH 3° “12M00O. + Co = (molybdenum blue) + Sn 


4) G 3] 


After 10 minutes of color developing its intensity 
was measured on spectrophotometer and calibration curve was 
then used to relate light transmission to actual concentra- 


tions. 


3. Test for Dissolved Oxygen. 


The Winkler method was used to determine dissolved 
oxygen levels in our samples. This method uses manganese 
sulphate (Mnso,, ) and the alkali-iodine reagent (NaOH + ee 
The oxygen present in sample oxidizes the manganese ion Mn 


to a higher valence and manganese oxide then precipitates 
as a brown flocculant substance. This part is so called 


fixation of oxygen. 
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The floc is allowed to settle and then concentrated sul- 
furic acid is added. Under resulting low pH conditions 


MnO. oxidizes I. to produce I,- 
MnO, + 217 + 4H’ = Mn” +1 


The iodine in the sample is then titrated with 0.025 N sodium 
thilosulfate? “Amount*of titrated iodine is’ equivalent to 
the amount of oxygen that was present in the water. 
4Na2550, a 21, = 2Na,S,0- + 4NaI 
Starches used aS uanelndiGateorecom thisatr eration «s.Added 


starch produces blue color which gradually disappears when 


alt gLodine wis titrated . 


4. Biochemical Oxygen Demand Test. 

Biochemical oxygen demand (BOD) is usually defined 
as the amount of oxygen required by the bacteria while 
stabilizing decomposable organic matter under aerobic 
conditions. The BOD test is widely used to determine the 


pollutional strength of domestic and industrial wastes in 
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86 
terms of the oxygen that they will require if discharged 
into natural stream. By its use it is possible to deter- 
mine the degree of pollution in stream at any Eimer 

The BOD test may be considered as a wet oxidation 
procedure in whieh the living organisms serve as the medium 
for oxidation of the organic matter to carbon dioxide and 
water. A number ofthe organisms responsible for the 
stabilization of organic matter in natural waters are 
forms native to soil. Theoretically an infinite time is 
required for complete biological oxidation of organic 
matter, but for all practical purposes, the reaction may 
be considered complete in 20 days. However, a 20 day 
period is too long to wait for results in most instances. 
It has been found that a reasonably large percentage of 
the total BOD is exerted in 5 days, consequently the test 
has been developed on the basis of a 5 day incubation period. 
It should be remembered, therefore, that 5 day BOD values 
represent only a portion of the total BOD, usually 70 to 
80%. 

The BOD test is based upon determinations of dis- 
solved oxygen levels. Samples whose BOD does not exceed 
7eitg / Ll ocenocercguite dilution, if they initially contain 


sufficient dissolved oxygen, The procedure for this 
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87 
type of sample is very simple. There are two or more 
samples taken for the test. The first is analyzed for 
dissolved oxygen immediately, the second is analyzed after 
incubation for 5 days at 20 degrees C. The 5 day BOD is 
then determined by subtraction of residual dissolved oxygen 


from the initial value. 


5. Laboratory Tests Results. 

All results of laboratory tests of water samples 
from the North Saskatchewan River taken at Vinca Bridge 
and Fort Saskatchewan Bridge are shown in Appendix C. The 
difference in chlorides in the river cross-section at Vinca is the 
best: indication: of the mixing occurring between Fort 
Saskatchewan and Vinca. The other tests as indicated by 
phosophates also confirm incomplete mixing, but because of 
their relatively low concentrations the uneven distribution 
throughout the river cross-section is far less apparent. 
Similar results are for the Total Carbon (TC), Dissolved 
Oxygen (DO), and Biochemical Oxygen Demand (BOD) tests. 

Laboratory test results from samples taken at Fort 
Saskatchewan Bridge show quite good mixing. There is no 
significant difference between the values of phosphates at 


Vinca and Fort Saskatchewan. Small difference in the DO 
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and BOD tests confirm DO utilization in the reach, which 
stays closed for the whole winter period with no reaera- 
HLone rossi ble. 

All laboratory test results are tabulated and to- 


gether with corresponding diagrams in the appendix, part C. 


= 


fw 
- 







bi 


dv <tisey mit (fb ip teeeeeeee ~ ine’ 


-s1p0e8. of dd behing sean okerne = 


. Pa : he Ve o& 
nop Bete stiiiet Say [TOU nay Ded leiesedat Lin” 
Fd ee oe 2 KE Teh. i en woynse nt oye i ae 


Ate? ps 





89 


TABLE Cl 


CHLORIDES CONCENTRATIONS AT FORT SASKATCHEWAN (mg/1) 





Distance wrom the might bank min feet 
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TABLE C2 


CHLORIDES CONCENTRATIONS AT VINCA (mg/l) 
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Distance from the right bank in feet 
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TABLE C3 


PHOSPHATES CONCENTRATIONS AT FORT SASKATCHEWAN (mg/1) 


a 





Distance from the right bank in feet 


Date ee oe ee 
100 200 300 400 
a 

03/10/71 O 0.36 0.28 0.22 ‘ate 
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P - polyphosphates 
T - total phosphates 
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TABLE C4 


CONCENTRATIONS OF PHOSPHATES AT VINCA (mg/1) 


er 


Distance from the right bank in feet 
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T - total phosphates 
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DAB oko 


DISSOLVED OXYGEN CONCENTRATIONS AT FORT SASKATCHEWAN AND 
VINCA CROSS-SECTIONS 















































el es Saskatchewan 
' Distance from the right bank in feet 
Date pata. 
100 200 300 400 
PEsL O77 1 oe 23 Leer e2e3 
ZOhIA fal IA 6) ra eee) 1230 
MG e/a iL ake fe) Deigees thalests! ive 
eA se ye rial a0 OF a) 9.8 Sieve: 
POF O37 bz O56 9.6 9.8 ome 
Tove Os) ae os Saskatchewan 
Date LOO 200 300 400 
eee 
08/03/72 522 3.4 B0 SU 
OZ O03 /e2 Ziel: Zee Dire al 20 
DEO. a) Vinca ee 
Date 100 200 300 400 550 
nS es ee ee ego 
7 LO sei i akib le? aba ae! fee 126.0 
26/2 fe). eae 1270 eee Lea I2Zne 
Wee Te eOrere LOe4 ith a ee ie 2 
08/03/72 9.0 Ss) 9.4 9.6 9.6 
O07 OS 7al 2 8.8 Oru Dee 9.6 SESS) 
Sa i elt tll a a A: a a 
BlO7D. =o avLica 
Date 100 200 300 400 50 
nner re 
08/0387 72 Bee Bele) AO Shs t ah 2 
ie) 7.08 Jai 2 Shige Sd) 26 2.4 22 









OO “4teRaAt 


GUA NAMI TALSAL : STAI imO WIS aiviosera 4, 
. eee 
e shad Pe a4 j elt * LETS | = ‘ — 

a - ao a ' i ot ne e be : 38:8 


im 27H SEAG : 
a tee oyna 





ad IT\Or\ er 
od [T\TL\be 
if IiTyer\al 
1c iT\e0\ee 
r St, co, ol 





oted 





—_ 


inf ev\eo,eo 
ET\EG\OL 





96 


‘ueMEYyo IeYsSeS JOY Je UOTIOOS-SsorH oy UT UoTINqTAISTd *q'O'd 


qeoy UT yuUeq JYSTAI oy. worZ DdURISTIG 








yueq 3391 OSS 004 00€ 002 OOT 
aaa oe aioe | ri fas a) | a aaa iat Ree 
ifs / Olen ee ee 
CL | Oy ee wee eee cae 
Saeed ~ 
a =e —— a 
SS See Se 


> €O anos 


a 


0 


ci 


uo TZeAQUDDUND 


qT /su 





oh 


CLC) Gee ey 


*"BOUTA JB UOTIDOG-ssoIg 9YI UT UWOTINGTIAISTq *a’oO’|g 


qaoy UT YUeq 4YSTA sy. WoAF soURrISTC 


yueq 3391 OTL 


eer StS era a ST ea | eS at aria nS RE A | re a age Aa a 4 ag 


007 O0€ 002 OOT 


AUG OO 


Sey 


Sf 


90 FNNDIA 


i) 
™N 


Oo 
(sa) 


fo) 


0°S 


i eS Sei eee eee ee Oe eee sa Oo 


UOTIEAIUIOUOY) 


T/ sui 





—_ 
x 


: = 
dned. ei Of' ozt ie we oes 
3667) @b tow 22 22G0) eeu) ssH0I4—7% 


aani¥..¢@ Lost ) «6 el a 














a vil | 

f x - 

- 7 ae oa 
ear | 


ce 
3 
alr 
- 








